Transformation from spots to waves in a model of actin pattern formation 
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Actin networks in certain single-celled organisms exhibit a complex pattern-forming dynamics 
that starts with the appearance of static spots of actin on the cell cortex. Spots soon become 
mobile, executing persistent random walks, and eventually give rise to traveling waves of actin. 
Here we describe a possible physical mechanism for this distinctive set of dynamic transformations, 
by equipping an excitable reaction-diffusion model with a field describing the spatial orientation 
of its chief constituent (which we consider to be actin). The interplay of anisotropic actin growth 
and spatial inhibition drives a transformation at fixed parameter values from static spots to moving 
spots to waves. 



Introduction. Dictyostelium discoideum {Dicty) is an 
amoeba known to generate spectacular patterns through 
organized multicellular aggregation fD. This organiza- 
tion is made possible by Dicty^s ability to move, which in 
turn is regulated by polymerization of the protein actin 
into oriented networks within individual amoebae [2^ . In- 
terestingly, these networks exhibit their own distinctive 
patterns. When treated with the drug latrunculin, actin 
networks within Dicty degrade, rendering the amoeba im- 
mobile. Upon removal of latrunculin, actin networks re- 
polymerize through a complex pattern-formation process 
that appears to consist of three stages O [4] , summarized 
in Fig.[l] First, immobile circular spots of actin form on 
the cell membrane. Second, spots acquire a persistent 
diffusive motion. Third, wave-like actin structures rem- 
iniscent of cell 'leading edges' appear and coexist with 
spots. As waves strike the cell periphery the amoeba 
recovers its ability to move. 

This paper describes a mechanism that might under- 
pin the remarkable dynamic transformation from static 
spots to moving spots to waves. Our starting point 
is the recognition that spots (both static and moving) 
and waves have been seen in chemical systems [5 , and 
that such patterns can be described by equations mod- 
eling reacting and diffusing chemicals [6]. Here we con- 
sider an excitable 'activator-inhibitor' reaction-diffusion 
model that describes, in different regions of its parame- 
ter space, stationary spots and moving waves. In order 
to interpret the chief constituent of this model (the 'ac- 
tivator') as substrate-bound actin in Dictyostelium, we 
equip it with an additional variable that describes local 
actin fiber orientation. We find that the resulting model 
exhibits a series of nonequilibrium transformations, at 
fixed parameter values, from stationary spots to moving 
spots to traveling waves. This transformation is driven 
by the interplay of inhibition (which permits localized 
spots) and the directional bias imparted to actin poly- 
merization by local fiber orientation, the latter emerging 
from a spontaneous breaking of fiber symmetry. Our re- 
sults support the idea that emergent actin patterns in 
vivo are fostered by an excitable medium [7^ 8 , and il- 
lustrate the dynamical richness accessible to an intrinsi- 
cally anisotropic chemical species that suffers inhibition. 



Models of inhibited but spatially isotropic chemicals can 
describe striking transitions between static and moving 
spots, but only as their parameters are varied [9]; tradi- 
tional models of anisotropic actin polymerization (tread- 
milling) do not possess as solutions static, size-limited, 
isolated spots. 



min 


13 min 


16min • 





FIG. 1: Dynamic actin structures in the substrate- attached 
cortex of a Dictyostelium cell, as visualized by TIRF mi- 
croscopy. LimEAcoil-GFP was used to selectively label F- 
actin. The cell was incubated for 15 min under 5 fiM latrun- 
culin A to depolymerize F-actin. The first frame (0 min) is 
taken 16 minutes after reducing latrunculin A concentration 
to 1 /iM; the first stationary actin spots are visible (scale bar 
10 /im). Spots become more numerous and mobile (13 min 
frame; selected traces over a period of 30 sec shown in red). 
At 16 min a prominent spiral wave appears (see dashed out- 
line). Bottom right: 4 frames of the counterclockwise-rotating 
spiral wave superimposed in different colors (red, green, blue, 
white; first and final images separated by 45 sec). 

We present our model below. We discuss the origin 
of spot stability, and demonstrate numerically that lo- 
cal fiber orientation can both induce a spot to move and 
distort a spot into a wave. The resulting dynamical evo- 
lution of the model captures several features of the spots 
and waves seen in vivo. 

Model. The FitzHugh-Nagumo equations [l0l[TT] stud- 
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ied by Vasiev [12] describe localized spots and waves in 
different regions of parameter space. We take these equa- 
tions as our starting point in an attempt to describe sim- 
ilar patterns seen in Dicty^ and augment these equations 
to account for the polarity of the fibers comprising poly- 
merized act in structures. We consider the evolution of 
scalar fields u{x^ t) and v{x^ t) on a two-dimensional sub- 
strate (cell membrane) according to the equations 

dtu{x,t) = V - J ^h{u) - pv ^ ^/2I\,r]{x,t)] (1) 
dtv{x,t) = S\/'^v^e{u-v). (2) 

Here J = \/u + Vqtu and h{u) = —ku{u — ui){u — uq). 
We shall set the parameters e, J, p, /c, uq and Ui by 
comparison with previous work [12 , and will explore the 
effect of varying Vq and T^. We consider u{x^t) to be 
proportional to the number of actin fibers per unit area 
(relative to a reference concentration) at substrate posi- 
tion X at time t. The field v{x^t) acts to degrade actin, 
and the vector field T{x^t) labels the orientation of the 
actin network. Equations ([T]) and ([2| (with Vb = = 0) 
describe classical excitable behavior wherein the auto- 
catalytic 'activator' field u and activator-suppressing 'in- 
hibitor' field V interact to generate localized patterns that 
may propagate spatially [13 . We consider the homoge- 
neous terms of Equation ([T]) to describe polymerization 
of actin at a rate proportional its local concentration, 
^poi = ^('^1 uq)u^^ and degredation of actin with rate 
^deg = —kuouiu. We regard the term —ku^ as a model 
of the concentration-limiting effect of steric hinderance. 
The field 7^ is a Gaussian white noise with zero mean and 
unit variance; Tu > quantifies the magnitude of this 
noise. 

We motivate the terms in Equation ([T]) describing 
spatial propagation of actin in the following way. We 
hypothesize (on the basis of bleaching experiments in 
Dicty [4 ) that stationary spots are composed of fibers ori- 
ented principally normal to the substrate (Fig.|2|a)). We 
assume that membrane-bound proteins (such as MyoB) 
recruit Arp2/3 complexes to the growing ends of fibers, 
and that these complexes in turn initiate the growth of 
daughter fibers at the spot periphery. We argue that 
such growth is isotropic, and the resulting propagation 
of material diffusive. By contrast, we expect that waves 
possess many fibers aligned in part parallel to the sub- 
strate and move chiefiy by treadmilling [2, J^. To de- 
scribe fiber orientation we have introduced a field T{x^t): 
vanishing r describes fibers pointing solely in the verti- 
cal direction, while nonzero r describes fibers with some 
component of orientation parallel to the substrate. We 
assume that lateral fiber orientation biases the direction 
of actin growth, which we model using the term Vo\/ -tu. 
Vq controls the rate of directed polymerization. We fur- 
ther hypothesize that orientation may be acquired spon- 
taneously, and therefore require that r evolve according 
to the equation 

r;^dtT{x,t) = --^ 1 d^xJ^r[r,u]^y/2T;ax,t). (3) 



Here J^r = h^r - ur'^ + r"^ + ai (V • rf + 0^2 (V x rf 
is a Landau-esque free energy density describing a field 
r that may order in the presence of the field u (pro- 
vided that the 'barrier' to initial ordering, the term in 
ho, is sufficiently small) and that has a tendency to align. 
Its degree of alignment is determined by a competition 
between the orientation-inducing terms coupled to ai,2 5 
and the orientation-destroying effect of the noise term 
coupled to T^. The field ^ is a Gaussian white noise 
with zero mean and unit variance. We set (arbitrarily) 
ai = a2 = 0.2, and regard T^- as our chief measure of 
network- alignment propensity. The kinetic prefactor F^- 
controls the rate of network alignment relative to that of 
network polymerization. 
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FIG. 2: Stationary and moving model spots, (a) Cartoon 
of the actin fiber structure we conjecture for stationary and 
moving spots in vivo, (b) Homogeneous effective 'free energy' 
density f{u) from our model (with uq = 0.05, k = 4.5 [12]) for 
values of actin polymerization rate ui below {ui — 1 solid 
line) at {ui — u\ ^ 1.06; dashed line) and above {ui — 1.1; 
dotted line) its first-order critical point; (c) nullclines of Equa- 
tions ([1]) and ([2]) under the same conditions, (d) ix-field pro- 
files {Id cuts through a 2d box of size 500^) at three differ- 
ent times ((ti,t2,t3) = (5,10,20) x 10^) from three differ- 
ent simulations (Vb = 0,0.012,0.06). We excite the medium 
at t = so that a spot nucleates, and we impose an ori- 
entation field T = x/2. The spot is stationary and stable 
for Vb = 0; the spot moves persistently without significant 
structural distortion at the two larger values of Vb- We show 
at ta the v-field profiles also (thin lines). Parameter values: 
Uq = 0.05, = 1, A; = 4.5, e = 0.5, 5 = 2.5, Tu = 0. 
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Equation ([2| describes the evolution of an actin- 
suppressing 'inhibitor field'. Actin recruits the agent of 
its destruction at rate eu; this agent degrades actin at a 
rate —pv^ and is itself degraded at a rate —ev. Possible 
biological origins for v include the concentration of actin 
severing proteins (e.g. cofilin) and actin capping pro- 
teins, or the state of hydrolysis of fibers. Here we simply 
regard 'i; as a coarse-grained agent of actin degredation. 
Our model assumes that actin polymerization is not in 
thermal equilibrium [2 . 

Intuitive explanation of spot stability. Other authors 
have demonstrated semi-numerically [12 and analyti- 
cally [151 [ISl [HI that equations similar to (1) and (2) 
with Vb = admit static spots as steady-state solu- 
tions in certain parameter regimes (a stationary spot 
profile is shown in Fig. [2] (d)). Here we follow ap- 
proaches detailed in Refs. [9l [161 (III [H] to put for- 
ward a simple physical argument for why such spots can 
exist. At steady state, and with coordinates rescaled 
such that (e/(^)V^, Equation ^ has solution 

v{x) = (27r)~-^/ dx'Ko{\x' — x\)u{x'). We expand this 
solution to second order in V and insert the resulting ex- 
pression in the steady-state version of Equation ([T]) . We 
regard the equation so obtained as the Euler-Lagrange 
equation of the 'free energy' functional 

{f{u)+g,{Vuf + c,{V^uf], (4) 



asu + a2{p)u ; = k/4:; as 



Fu[u] = j (fx 

where f{u) = a/^u^ 
k{uo + Ui)/3; a2{p) = {kuiUo + p)/2; gp = {e/d - p) /2; 
and Cp = p/2. The homogeneous component f{u) de- 
scribes a first-order phase transition from an empty sub- 
strate to an actin-covered substrate as the parameter ui 
is increased {ui is roughly proportional to the mean-field 
actin polymerization rate); bulk phase coexistence occurs 

1/2 

when m = Ui = ^uq -|- | (i^q + 8p/k) . For the param- 
eters of reference [12 (p = uq = 0.05 and k = 4.5), we 
have Ui = 1.0632; in that paper, ui was set to unity, im- 
plying that the system considered there was numerically 
close to coexistence. We show f{u) (and its correspond- 
ing nullclines) in this regime in Fig.|2|b) (Fig.|2|c)). Note 
that the activator-inhibitor coupling p infiuences the lo- 
cation of the mean-field critical point, but not the nature 
of the transition. 

The space-dependent terms of Equation Q admit 
modulated patterns when p > and e/5 < p; in this 
paper we perform simulations with p = l,e = 0.5 and 
6 = 2.5. It is instructive to recognize that while versions 
of Equation Q appear in many different settings [17 , 
the inhibitor-induced spatial modulation of the u field 
closely resembles the physical mechanism by which sur- 
factant induces microphase separation of oil in water [TF . 
One might therefore consider the stable spots of the 
FitzHugh-Nagumo equations studied by Vasiev [12 to be 
akin to drops of surfactant-coated oil in water. Moreover, 
extremizing Equation Q in circular geometry for an as- 
sumed density profile (adapting the approach of Ref. p[8] ) 
reveals a characteristic spot size. While a local approxi- 



mation such as Equation Q possesses limited predictive 
power [T7| , it does capture the correct trend of variation 
of spot radius with model parameters (such as e), and 
provides an intuitive explanation for why static spots of 
well-defined size exist in this system. 
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FIG. 3: Model spot-to- wave transformation, (a) An initially 
static spot moves under the influence of a fixed orientation 
field T = x/2 (Vb = 0.24) and deforms into a structure re- 
sembling a traveling wave (profile and image shown 1.05 x 10^ 
time units after application of the field) . (b) We quantify this 
distortion mechanism by plotting the radius of gyration Rg 
as a function of time for spots exposed to similar fixed r and 
various values of Vb- (c) Time-ordered snapshots of i^(x) from 
a simulation with fluctuating r starting from a single spot in 
the center of a box. The spot acquires orientation sponta- 
neously (spot arrow labels spot- averaged r), executes a per- 
sistent random walk (numbers identify periodic box crossings; 
trajectories and arrows record spot motion), and begins to de- 
forms into a wave at t ~ 24 x 10^. (d) A similar simulation 
with nine initial spots mimics the transformation from static 
spots to moving spots to waves shown in Fig. [l] 

Effect of fiber orientation upon spots. To determine the 
effect of T upon the behavior of spots, we performed nu- 
merical simulations of our model for a parameter set that 
admits, at steady state, immobile spot solutions when 
Vo = {uo = 0.05, ui = l,k = 4.5, e = 0.5, 6 = 2.5 [12]). 
We integrated Equations (1-3) using the second-order 
Heun algorithm [19 on a square grid of lattice constant 
0.4 with timestep 0.01 [12 . We used a 9-point Lapla- 
cian [20]. In Fig. |2|d) we show that when r and Vq 
are fixed, a spot moves persistently. Moving spots in our 
model behave similarly to the moving spots described in 
Ref. [9 , colliding both inelastically (at high speeds) and 
elastically (at low speeds). When Vq is made sufficiently 
large, a spot deforms into a localized structure (Fig. 
|3ja,b)) resembling the traveling waves seen in the regu- 
lar FitzHugh-Nagumo model at parameter values distinct 
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from those at which spots are found [12]. This exploita- 
tion of a spot's apparent underlying instability to wave 
deformation, which embodies the notion of nonequilib- 
rium pattern control discussed in Ref. [21 , requires per- 
sistence of orientation. In the full model, when r evolves 
according to Equation (|3| (with noise), the transforma- 
tions from static spot to moving spot and from moving 
spot to wave can occur on a broad range of timescales. 
In Fig. |3|c) we show configurations of u{x) from a simu- 
lation (r^ = l,ho = 0.225, Tu = 0,Tr = 0.015, = 0.1) 
that starts from a single spot at the center of a box of 
size 450^, periodically replicated in imitation of bulk sur- 
roundings. The spot spontaneously acquires orientation 
(at t 1.3 X 10^), executes a persistent random walk, 
and eventually deforms into a wave. Spot trajectories 
are shown in the second and third panels, with periodic 
boundary crossings labeled in the order they occur. In 
Fig. |3|d) we show images from a similar simulation with 
nine initial spots. A phase of autonomous mobile spots 
is supplanted at later times by an organized collection of 
waves (moving from left to right), mimicking the trans- 
formation seen in Dicty. Indeed, the model's dynamics 
captures several qualitative features seen in vivo [l] 0]: 
spots nucleate spontaneously (when > 0.05); spots 
spontaneously (and autonomously) acquire orientation; 
oriented spots execute persistent random walks; spots 
split (cf. [22^) and merge; moving spots deform into 



larger, wave- like structures on a broad range of deforma- 
tion times; waves and spots may coexist; and large struc- 
tures degrade or limit the growth of neighboring smaller 
structures. 

We have shown that a simple and physically-motivated 
modification of the FitzHugh-Nagumo equations can 
mimic the interconversion and coexistence of distinct dy- 
namical motifs seen in Dicty ostelium. Spatially local- 
ized patterns in reaction-diffusion systems have recently 
received considerable attention [6], and the model we 
present demonstrates the unusually rich dynamics acces- 
sible to an intrinsically anisotropic chemical species. Our 
work suggests that certain cytoskeletal dynamics can in- 
deed be caricatured by simple reaction-diffusion models 
embodying the idea of excitability [71 [8], supporting a 
physical picture (that one might dub the excito skeleton) 
in which cytoskeletal pattern formation is driven by the 
self-organization of excitable solitons of actin and its at- 
tendant proteins, rather than being orchestrated solely 
by biophysical signaling pathways. 
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